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Abstract 

We study the four-body decays J/ip — > NN meson meson using a chiral unitary ap- 
proach to account for the meson meson final state interaction (FSI). The calculation 
of the J /ip —>■ NNtc + tt~ process properly reproduces the experimental data after tak- 
ing the FSI of mesons and the contribution of intermediate p meson into account. 
The isoscalar resonances a, / (980) and the isovector resonance a (980) are generated 
through the FSI of the mesons in the channels J/ip —>■ NNtt°tc and J/ip — > NN7r°r), 
respectively. We also calculate the two mesons invariant mass distribution and the 
partial decay width of J/ip — > NNK + K~ and J/ip — > NNK°K°, on which there is 
still no experimental data available. 



In past decades there have been many efforts in understanding the nature of hadronic 
interactions. Such efforts follow basically two strategies. One of them describes the 
hadronic interactions in terms of the quark-gluon structure of the hadrons, the un- 
derlying theory of which is QCD. The other one, chiral perturbation theory (ChPT), 
deals directly with mesons and baryons at low energies. It incorporates the basic sym- 
metries of QCD into an effective Lagrangian expanded in powers of the momentum 
of the hadrons with which one performs the standard field theoretical calculations for 
meson-meson [1] and meson-baryon [2, 3, 4] interaction at the lowest orders. ChPT 
provides an elegant, systematic and technically simpler way to make predictions in 
hadronic processes at low energies, where QCD becomes technically unaffordable, and 
has achieved many successes. The obvious drawback of ChPT is its limited range of 
convergence. For instance, for meson meson interaction this limitation appears around 
500 MeV where the a pole shows up. Plain ChPT can do little for the investigation of 
the interesting resonances that occur in meson spectroscopy. 

Recently, a chiral unitary coupled channels approach, which makes use of the stan- 
dard ChPT Lagrangian's together with a expansion of Re T -1 instead of the T matrix, 
has proved to be very successful in describing the meson meson and meson baryon 
interactions in all channels up to energies around 1.2 GeV in meson meson and 1.6 
GeV in meson baryon interactions [5, 6, 7]. With the coupled channels Bethe-Salpeter 
equation, a technically simple calculation was done in [5] for the S-wave meson meson 
scattering using the lowest order ChPT Lagrangian as the source of the potential. It 
was shown there that the effect of the second order Lagrangian can be appropriately 
incorporated and the resonances a, / (980) and a (980) were generated as poles of the 
S-wave amplitudes. This approach has been used in the investigation of the final state 
interaction (FSI) of mesons in several decay processes in order to get a better under- 
standing of the nature of meson resonances. In this aspect the authors of [8, 9, 10] 
studied the FSI in radiative decays, getting a very distinct peak for the resonance 
/o(980) in the titt invariant mass distribution of the <fi — > 7r°7r°7 decay and a dominance 
of the a (980) for -> 7^7. 

In this work we investigate the four body decay of J/ip into NN meson meson. 
The J/ip decay with NN in the final state has already attracted some theoretical 
attention, albeit with only one meson in the final state [11]. We treat the mesons FSI 
with the techniques of the chiral unitary approach. The final meson pairs we consider 
here are it + it~, 7r°7r°, K + K~, K°K° and Tr°r] which can be expected to produce meson 
resonances through FSI. The J ftp has a mass of 3.097 GeV, but in the decays studied 
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here a large fraction of this energy is lost creating the NN pair. The remaining energy 
for the pair of mesons falls well within the range of the energies where the chiral unitary 
approach has been successfully applied. 

In Fig. 1 we present the diagrammatic description for the J/ip — > NNMM decays 
including the meson meson FSI. 

In order to construct the amplitudes for J/ip — > NNMM we take into account 
that J/ip can be considered as a SU(3) singlet [12]. Then we write the most general 
NNMM Lagrangian of SU(3) scalar nature without derivatives in the fields. We have 
the following possible structures for the Lagrangian: 

d = g Tr[B'fB&&\il li , C 2 = g Tr[B'f&B§\V 
C 3 = g TriB-fQQB^p, C 4 = g Tr[B^B]Tr[^]^^ (1) 

with $, B the ordinary SU(3) matrices for pseudoscalar mesons and | + baryons, re- 
spectively, the J/fy field and g a constant to provide the right dimensions. We then 
take the Lagrangian of our problem as a linear combination of C a , a — 1, 2.. 4, 

C = "£x a C a . (2) 

o=l 

This leads to the vertex 

Vi = -d g u(p')Yv(p)e,(J/^, (3) 

where we have already specified that we have a baryon anti-baryon production, rather 
than the baryon destruction and creation that one has for the meson baryon amplitude. 
We shall only study the cases where we have pp and nn production. Thus we can 
have J/ip — > NNii + 7r~, 7r 7r°, K + K~, K°K°, 7i°r] and we will call these channels 
2=1, 2.. 5, respectively. The Cj coefficients are given in Table I. In the rest frame of 
J/ip, 6r(J/ip)=0 for the three polarization vectors and eq.(3) in the non-relativistic 
approximation for the nucleons can be written as 

V i = C i gff-^J/ 1 p) (4) 

with a the standard Pauli matrices. 

Table I : Cj coefficients in eq.(4) 





ppn + n 


pp7T 7l 


ppK + K~ 


ppK°K° 


ppn !] 


Ci 


x 3 + 2x 4 


x 3 + 2x 4 


x\ + x 3 + 2x 4 


x 4 + 2x 4 


^(-x 2 + x 3 ) 






nn7i°n 


nnK + K- 


nnK°K° 


nnn°i] 


Ci 


x 3 + 2x 4 


x 3 + 2x 4 


Xi + 2x 4 


xi + x 3 + 2x 4 


^(x 2 -x 3 ) 
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Note that the 7 M in the Lagrangian is of order 1 for NN production while it is of 
order q/M B for J/vpB — > BMM. The structure of eq. (2) is the only one without 
momentum dependence and in the present problem, where there is relatively small 
phase space for the final particles, one might think that it accounts for the largest 
part of an amplitude admittedly more complicated than the one of eq. (2). The fact 
that one has unknown parameters in the theory gives it a flexibility to approximate 
realistic amplitudes when fitting the parameters to experimental data. This approach 
is also used in [13] in the J/ip 07r°7r° decay where a simple structure of the type 
€ fi( < f>) e,l (J/i ; ) is adopted without derivatives. Some arguments are given there on why 
other structures would produce a weak s dependence and could be reabsorbed in this 
simple structure. Similar arguments could in principle be advocated here. 

The FSI in Fig. 1 is given by the meson-meson amplitude originated from the 
lowest order chiral Lagrangian with the chiral unitary approach [5], with which the 
amplitudes of decay channel % can be written as 



corresponding to the loops with the intermediate propagator of the i — th meson pair in 
which P is the total four-momentum of the meson-meson system with masses and 
m i2 and q is the four- momentum of one of the intermediate mesons, the loop integration 
variable which is regularized with a cut off \q\ < q max and q max =1030 MeV [5]. The 
matrix tij, which is actually a 5 x 5 symmetrical matrix, accounts for the meson-meson 
amplitudes between the i — th and the j — th meson pairs. 

We consider first the meson meson S-wave interaction and we only need the meson 
meson S-wave amplitudes in the loops of Fig. 1. By using the isospin Clebsch-Gordan 
coefficients the physical amplitudes Uj needed here can be written in terms of the 
isospin amplitudes of [5] : 



U = ^ + E Wit 



(5) 
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with G a diagonal matrix, with matrix elements 






ti4 — tn+n-^KOR — —f^tinT^KK 



■S3 

2 ,7=0 



^22 — ^TrOTi-'W'Tr — ~^^tttt — >7V7T 

1 



hs — tn°TT°^K+K ~~7kKtt^KR 



t 2 4 — tnOnO^KORO ~ ~7^K^KK 



1=0 

^33 = tK+K-^K+K- = T^KR^KK + ^KR^Kr) 
t M = t K+K -_> K OKO = — {tfcR^KK ~ ^KR^KR) 
^35 = t K +K-^rt = -^^KR^ttti 
tu = t K ORO^ K ORO = -^kR^kR + ^K~R—*Kr) 
t45 = t K ORO^O v = -^^KR^-q > 

^55 — t-n-O^^Or! = t^^^, (7) 

in which the small 1=2 amplitudes are neglected. In eq. (7), as in ref. [5], we use 
the phase convention \n + >= — 111 > and \K~ >= — \\ — | > in isospin states and 
a factor \[2 for each 7v + n~ and 7r°7r° state is introduced to compensate the unitary 
normalization of the isospin amplitudes which includes an extra factor 1 / \[2 for each 
of the isospin ixix states. 

There is an important point worth mentioning concerning the off shell part of the 
meson meson amplitudes in the meson loops. In eq. (5), the meson meson amplitudes 
in the loops are taken on shell. It was shown in [5] that the contribution of the off shell 
parts could be reabsorbed into a redefinition of coupling constants in meson meson 
scattering. Analogously, for the vertex V, h which has a different structure than the 
meson meson amplitude, the contribution of the off shell part in the first loop could be 
absorbed by renormalizing the coupling constant g in eq. (1) since the integration for 
the off shell part in just one loop has the same structure as the tree diagram [5]. 

So far we have dealt with mesons in S-wave. We now turn to the P-wave. There 
is no 7r°7r° state in P-wave and the n°r] system also couples extremely weakly to P- 
wave [16], but the n + n~ state can be in P-wave driven by the p meson. Due to 
the strong coupling p° — > 7r + 7r~, the contribution of pn + n^ coupling to the decay 
J — > NNir + ir~ depicted in Fig. 2 should be considered. Analogously to the ppNN 



4 



coupling discussed in [17, 15], the J/ipNNp vertex has the structure 



-Hj p nn = ~i 9j p nn e \p) ■ KJ/^V , ( 8 ) 
where the index i stands for isospin, and the pirn coupling has the structure 

-it^ = i r ^e lx {p){p- P r , (9) 

in which p p and p'^ are the four- momentum of n + and 7r~ respectively. We take the p 
mass m p = 770 MeV and G v = 56 MeV. It can be seen from eqs.(4) and (8, 9) that 
there is no interference between the amplitudes described in Fig. 1 and Fig. 2, and 
hence we have the total strength for the decay J/ip — > NNn + n^ by simply summing 
the strengths of the mechanisms of Fig. 1 and Fig. 2 with the intermediate p° 

EE \Ti\ 2 = EE M 2 + EE K\ 2 (10) 

with 



2 1 2 ,m n G\ 



EE l*!pl ~~ ^9j p nn( 
where 



P 



Ml — m 2 + iY p (Mi)mr 



2 

2 /l™2 



(Mj*-4m£) (11) 



W)= r ^g» P=(j-^. Pb = (^-"*)* (12) 



P ^,2' 



with Mi being the 7T7T invariant mass and V p = 149.2 MeV the experimental p width. 

We define the amplitudes Tj=£j for i=2..5, which implies that we neglect the p 
term in tt°tt°, K + K~~, K°K° and ir°i] production. The p does not couple to 7r°7r° for 
symmetry reasons. Its coupling to pairs of mesons can be obtained from the chiral 
Lagrangian of [18], where we find that it does not couple to ir°r] but it couples K + K~ , 
K°K° with strength ±, -\ that of 7r + 7r . We have evaluated the contribution of 
the p term to K + K~ and K°K° production and find it negligible compared to the 
contribution coming from the tj amplitude. 

Considering that Tj is a function of the single invariant masses of the final mesons, 
the meson meson mass distribution of the four body decays can be reduced to 



2 ™2 ™2 _ 

' • (13) 
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where the Heaviside step function 6{x) and the Kaellen function X(x, y,z) have been 
used, with A given by 

A = ^l M J + 2M * + 2E (pW) - M i - 2Mj(E(p) + E[p'))\ , (14) 

where Mj is the meson meson invariant mass of the i-th decay channel, Mj and M are 
the J/ip and nucleon masses, respectively, and E(p) and E(p') are the energies of the 
final nucleons with 3-momentum p and p', respectively. 

It can be seen from table I that the Lagrangian C 2 in eq. (1) does not contribute 
to the nir spectrum. We define the ratios r\ = with i=l, 2, 4. The strength and 
the shape of the nn spectrum at lower energies are mainly determined by the constant 
9a = (2^3 + 2x4)g, while the shape at higher energies is adjusted by the ratio r 4 . It 
is worth noting that the ratio n also influences the shape of the nn spectrum at 
higher energies, but its contribution could be included in the variation of r 4 to fit the 
experimental data. The parameter r 2 does not change the nn spectrum but plays an 
important role for the n°r] and KK spectrum. 

Fitting the data, we find strong restrictions for the r 4 parameter. Indeed, we can 
reproduce the data only for values of r 4 in the range —0.3 < r 4 < 0.3. For values 
of r 4 outside this range a large signal for the / (980) excitation, incompatible with 
experiment, is obtained. A best fit to the data is obtained with r 4 = 0.2, in which no 
peak for f excitation is seen, and with r = —0.27, where a signal for f excitation 
compatible with experiment, assuming the single datum in the / (980) peak is not a 
statistical fluctuation, is obtained as one can see in Fig.3. The values for the coupling 
constants are g a — 1.1 x 10~ 6 MeV~ 2 and g JpN fj = 7.1 x 10~ 5 MeV~ 1 for the case of 
r = 0.2 and g a = 1.2 x 10~ 6 MeV~ 2 for r = —0.27 and the same value for g JpN N- 

By using the above values of coupling constants g a and g.j P NN, we derive the con- 
tribution to the width from the mechanisms in Fig. 1 and p° to be 4.4 x 10~ 4 MeV and 
9.1 x 10~ 5 MeV, respectively. Then the width, which is the summation of the two contri- 
butions, yields 5.3 x 10~ 4 MeV, or equivalently a branching ratio 6.0 x 10~ 3 , consistent 
with the experimental branching ratio (6.0 ± 0.5) x 10~ 3 of the decay J ftp — > ppir + ir~ 
[19]. The p branching ratio that we obtain is 1.05 x 10 -3 . This seems to be in conflict 
with the present data in the particle data group which quotes a value smaller than 
3.1 x 10~ 4 . However, one should note that an ordinary fit to the 7r + 7r~ data with a 
background following phase space and a p contribution would lead to a smaller p con- 
tribution than the one we obtained here. Yet, as we discuss below, the contribution to 
7T + 7r~ of the chiral terms does not follow phase space, and this is one of the interesting 
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findings of the present work. Indeed, in Fig. 4 we split the contribution to the 7r + 7r~ 
production from the mechanisms of Fig. 1 (for r4=0.2) into the Born contribution 
(Fig. la) and the meson rescattering terms (Fig. lb, c.) plus the interference of 
the two terms (dotted, dashed and dash-dotted curves), respectively. The contribution 
of the Born term in Fig. 4 follows phase space, which, as we can see in the figure, 
is quite different from the final invariant mass distribution. The contribution of the 
rescattering terms, which involves the meson meson interaction in S-wave, accounts 
for a production seen as a broad bump, which does not exactly follow phase space, 
and a small peak that accounts for / (980) excitation. The sum of these two contribu- 
tions would not differ very much from phase space, apart from the small /o(980) peak. 
However, the interference term is peculiar and changes sign around 500 MeV and the 
resulting invariant mass distribution is much narrower than the phase space or the a 
shape and loses strength in the region of the p excitation with respect to the phase 
space distribution. The strength around this region is filled up by p meson production 
in our approach and this leads to a bigger contribution than the nominal one in the 
PDG. 

By using the parameters which reproduce the experimental data for J/ip — > ppn + 7i~ 
decay, we can give a prediction for J/ip — > ppn 7i° decay. In the model employed here, 
one can see from table I and eq. (7) that the Born term for J/ip — > ppir°ir° is the 
same as that of J/ip — > ppir + iT~, and the rescattering terms are also the same. This 
means that the the shape of the 7r + 7r~ spectrum and the 7r°7r° spectrum should be 
the same, apart from the p contribution to 7r + 7r~. We present the n°TT° spectrum for 
J/t/j — > ppTc°Tc° in Fig. 5 (the strength is divided by a factor of 2 because of the identity 
of the two 7r°s). The width of the decay J/ip — > pp-K°-K° comes out to be 2.2 x 10~ 4 
MeV, which corresponds to a branching ratio of 2.5 x 10~ 3 . It would be interesting to 
measure this decay to search for the different strength of the distribution compared to 
the decay J/ip — > ppn + it~ . In our model, the a resonance has a clearer manifestation in 
these reactions. Indeed, the peculiar shape of the interference term between the V and 
the Born term for J/ip — > ppix+ix^ as shown in Fig. 4 has more drastic consequences, 
narrowing the tttt distribution and shifting the peak to small invariant masses, which 
are clearly visible in the data, although results with better statistics would permit to 
further check this observation. 

Our next step is the investigation of the J/ip — > pp7r°r] decay. According to table 
I, this channel is determined by the parameters r 4 and r 2 . Since r 2 does not influence 
the behavior of the n + 7i~ spectrum, which provides the only experimental data to fix 
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the parameters in our model, it is a free parameter which we vary within a reasonable 
range to see its influence on the n r] spectrum. Our observation is that there is always 
a peak that accounts for the a (980) meson in the 7i r] spectrum, while its branching 
ratio varies within a wide range from 1CT 4 to 1CT 3 for a range of r 2 between -1 to I. As 
an example, we take r 2 =0 to give a qualitative impression on the ir°r} spectrum in Fig. 
5. Since the a (980) stands firmly in our model, it would be interesting to search for 
it in the n r] mass spectra of this decay, which certainly will be constructive in fixing 
the parameters in our model. 

We also evaluate the decays J/ip — > ppK + K~ and J/ip — > ppK°K°, which also 
depend on the parameters r 4 and r 2 . As for n°rj, we also take r 2 =0 to qualitatively show 
the KK spectrum in Fig. 5 (the K°K° spectrum for r 4 = 0.2 and the K + K~ spectrum 
for r 4 = —0.27 are not shown because they are very small). It can be seen that the KK 
strength depends much on the model parameters while its shape is rather independent 
of them. By changing the model parameters we see that the branching ratios for these 
decays remain quite small compared to the other decay channels considered here. Since 
there is not yet experimental data available for the two decays, a future experimental 
measurement of these channels would also be most welcome. 

The method used here allows one to relate the J/ip — > nnMM with J/ip — > ppMM. 
We have also evaluated mass distributions for nfui+n^ and find the shape and strength 
very similar to the one of pp7v + n^ . This would provide an extra experimental test of 
the present theoretical approach. The branching ratio for J/ip — > nnn + 7r~ comes out 
to be 6.0 x 10" 3 , which is consistent with the PDG data (4 ± 4) x 10~ 3 [19]. 

In summary, within reasonable hypothesis, we have obtained a structure for the 
amplitudes of J/ip — * NNMM process, with still sufficient freedom, which allows 
us to get a good fit to the ir + n~ spectrum. The data does not fully determine the 
parameters of the model but puts some constraints on them which allow us to evaluate 
the uncertainties in the predictions. These show up as a broad range of values for the 
-k°t] branching ratio and KK channel. Yet, issues like the interference of the tree level 
and the V meson part of the amplitude, or the p contribution, remain unchanged 
within the freedom of the model and hence stand on firmer grounds. Experimental 
measurements of channels discussed here are certainly necessary for further progress 
in the understanding of the dynamics of the processes. 
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Figure 1: Diagrams for J/ip — > NNMM decays including the meson meson final state 
interaction. 

7l + 




Figure 2: The diagram for J/ij; — > NNtt + tt decay with a intermediate p meson. 
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Figure 3: (Color on line) The tttt invariant mass distribution for J/ip — > ppir + ir~ decay com- 
pared to experimental data with = 0.2 and -0.27. The solid line is the sum of contributions 
from the chiral terms and intermediate p meson, the dashed and dotted lines are the results 
of chiral terms and p contributions, respectively. Experimental data are taken from [20] 
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1 100 



M, (MeV) 

Figure 4: (Color on line) Decomposition of the contribution of the mechanisms of Fig. 1 
to the 7T + 7T™ spectrum for J ftp — ► ppir + Tr~ decay with = 0.2. The solid line is the total 
result from the chiral terms. The contribution from the Born term, rescattering terms and 
the interference between the two terms are plotted in dotted, dashed and dash-dotted lines, 
respectively. 
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